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Abstract 


In part I of the present study, a two-dimensional isothermal mechanistic model of a composite yttria/scandia stabilized zirconia anode-supported 
multiple layers solid oxide fuel cell (Ni-YSZ|Ni-ScSZ|ScSZ|LSM-ScSZ) was developed, calibrated and validated using a button cell setup. In 
this paper, the model was further validated for extended operating conditions of fuel and oxidant compositions, which confirmed the robustness 
of the model. Since an important attribute of this model is its ability to accurately describe the effect of cell microstructure on cell performance, a 
large section of part II is devoted to sensitivity analysis of various microstructure parameters (cathode area, porosity, layer thickness and particle 
radius) on cell performance. In particular, the simulation results have shown how the button cell microstructure should be modified to optimize 
its performance. The proposed model could be used as a bridge between micro-scale electrode models and macro-scale unit cell simulation, thus 


being a useful tool for electrode design. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


In part I of this study, a review of SOFC modeling 
approaches has been presented [1] and it was shown that high 
SOFC performance relies on superior electrochemical reac- 
tions and mass transport processes. Since experimental studies 
on SOFC are expensive, time-consuming and labor-intensive, 
quantitative mechanistic models for the cell PEN (positive 
electrode|electrolyte|negative electrode) structure are essential 
for improving SOFC designs. In part I, a two-dimensional 
isothermal mechanistic model of an anode-supported Ni- 
YSZ|Ni-ScSZ|ScSZ|LSM-ScSZ multiple layers solid oxide fuel 
was developed to describe the intricate interdependency among 
the ionic conduction, electronic conduction, multi-component 
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species transport, electrochemical reaction processes and elec- 
trode microstructure. The model is based on the following 
assumptions: (1) the PEN operates at steady state; (2) all gas 
mixtures are considered as ideal gases; (3) the reaction active 
sites are assumed to be uniformly distributed in each electrode 
layer. The two conducting phases (electronic and ionic) are 
considered to be continuous and homogeneous in each elec- 
trode layer; (4) the temperature is uniform in the PEN, thus 
the model is assumed isothermal and all physical properties 
are evaluated at a fixed cell temperature; (5) convection flux 
is neglected in the porous electrode compared to diffusion. 
Pressure gradients in the porous electrode are also neglected; 
(6) the boundary conditions for potential and species con- 
centrations are assumed uniform at the electrode/gas channel 
interface. 

Due to symmetry, a two-dimensional axial symmetrical 
coordinate is adopted. Refer to Refs. [1,2] for a detailed 
description of the governing equations. The calculations were 
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performed using the finite element commercial software COM- 
SOL MULTIPHSICS®, Version 3.2. 

The model was calibrated and validated by comparing var- 
ious overpotentials (anodic, cathodic and ohmic) determined 
experimentally with those determined from simulation. 

Beside the development of a comprehensive validated model, 
another objective of this work was to study the effects of various 
operating and microstructure parameters on the cell perfor- 
mance. Consequently, in the present paper, the model was further 
validated for extended operating conditions (different fuel and 
oxidant compositions). 

In addition, the present model is a continuum micro-scale 
model [3,4] which can, not only predict the detailed marco-scale 
cell performance, but also can predict the effects of micro- 
structural electrode parameters, such as electrode thickness and 
porosity. Thus, this model serves as a bridge between micro- 
scale electrode models and macro-scale unit cell simulation, 
which is an important attribute of this model as a design tool. 
Model parameters such as cell structure parameters, chemi- 
cal kinetics parameters, and transport properties were derived 
from experiments or published data. Note that for all simula- 
tion results presented here, all model parameters used are those 
that were determined in Ref. [1]; no additional tuning were 
performed in part II. 


2. Model validation and discussion for extended 
operating conditions 


2.1. Experimental and operating conditions 


The test setup and characterization techniques were presented 
in part I of the present study [1]. In Ref. [1] only one binary 
mixture of H2—-H20 was considered (4.2% H20 from humidifi- 
cation at room temperature). This binary mixture is referred to 
in this paper as the base case (humidified H) with air at 1 atm 
for temperatures between 750 and 850 °C). In the present study, 
operating conditions and validation of the model are extended to 
various water contents and to ternary mixtures where H2—H20O is 
diluted with N2. In addition, oxygen, instead of air, was also con- 
sidered as oxidant. The flow rates of fuel and oxidant were kept at 
50 sccm for all experiments. The polarization curves were mea- 
sured using four-probe method. Each operating variables was 
varied as follows: 


e cell temperatures: 750, 800 and 850°C; 

e oxidant: oxygen or air; 

e water content in the fuel: the fuel gas was saturated at 30°C 
(4.2% H20), 40°C (7.3% H20) and 50°C (12.4% H20); 

e Nọ dilution: N2 content in the dry fuel mixture: 0%, 20%, 
40%, 60% and 80% 


2.2. Effect of oxidant composition 


The SOFC performance can be improved by using pure O2 
instead of air as oxidant. Fig. 1 depicts the modeling and exper- 
imental results at different temperatures when pure Op is used 
as oxidant. Note that the curves are shown for current densities 


1.2 


1.0 


0.8 


0.6 


Cell voltage/V 


0.4 


@ Base case, experimental (Humidified H2/Air, 800°C) 
other I-V curves: Humidified H2/O2 


0.2 o 850°C, Experimental; —s— 850°C Simulated 
A 800°C, Experimental; —4— 800°C, Simulated 
* 750°C, Experimental; —*— 750°C, Simulated 


0.0 


0 2000 4000 6000 8000 10000 


Current density/((Am”) 


Fig. 1. Modeling and experimental polarization curves when using O2 as oxidant 
at 750, 800 and 850°C. Constant H20 concentration of 4.2%. 


up to 6500 A m7? because 1 A is the maximum current allowed 
with the electrochemical workstation (IM6e, Zahner-Elektrik 
GmbH, Germany), which corresponds to a current density of 
around 6500 A m~? for a cathode area of 1.54 cm’. 

It can be observed from Fig. | that the modeling results agree 
well with the experimental data at 750 and 800 °C, but somewhat 
underestimate the cell performance at 850°C. 

For comparison, the cell performance for the air case at 
800°C is also shown in Fig. 1. It is seen that the cell voltage 
at 800 °C is higher with pure O2 than with air. At a current den- 
sity of 6000 A m7”, the difference in cell voltage is about 80 mV. 
Fig. 1 also shows that improvement in cell performance due to 
the use of pure oxygen increases as the current density increases. 


2.3. Effects of fuel composition (mixtures of H2/H2O/N2) 


2.3.1. Effects of water content 

Fig. 2 shows experimental and simulated cell performance 
data using air as oxidant when the fuel gas was humidified at 
30, 40 and 50°C, corresponding to a water content of 4.2%, 
7.3% and 12.4%, respectively. Both experimental and simulated 
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Fig. 2. Effects of water content on the cell performance. 
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polarization curves show that, as the water content increases, 
the cell performance deteriorates. The modeling results agree 
well with the experimental data when the fuel gas contains up 
to 7.3% water. For higher water content (e.g. 12.4%), the mod- 
eling results overestimate the cell performance, especially at 
higher current density. As discussed in Refs. [2,5], at higher 
current densities and higher product water concentrations, the 
competitive adsorption and surface diffusion mechanism may be 
responsible to higher interface diffusion resistance. However, the 
current model does not take into account the effects of compet- 
itive adsorption and surface diffusion, which might explain the 
difference between experimental and simulated results at higher 
water content and higher current densities. 


2.3.2. Effects of nitrogen content 

Fig. 3 shows the experimental and simulated polarization 
curves for various mixtures of H2/H2O/N2 at 850 and 800 °C. In 
all cases, the water content was kept constant at 4.2%. Nitrogen 
does not contribute to the electrochemical reactions. However, 
not only N? could affect gas transport processes, but also, in the 
experiments shown in Fig. 3 the H20 concentration was constant 
(constant humidification temperature) and therefore the H2/H2O 
ratio was changed as the N2 concentration varied. It should be 


Temperature: 850°C 
Oxidant: Air 


Cell voltage/V 


95.8% H2,Experimental 

76.6% H2,19.2% N2,Experimental 
57.5% H2,38.3% N2,Experimental 
38.3% H2,57.5% N2,Experimental 
19.2% H2,76.6% Nz,Experimental 
Simulated 


“0 2000 4000 6000 8000 10000 


Current density/(Am”) 


95.8% H2,Experimental 
z 76.6% H2,19.2% N2,Experimental 
1.08 57.5% H2,38.3% N2,Experimental 
ý S 38.3% H2,57.5% N2,Experimental 
19.2% H2,76.6% N2,Experimental 
0.8 ẸŞ >S N g Simulated 


0.6 


Cell voltage/V 


0.4 
Temperature: 800°C 
Oxidant: Air 


0 2000 4000 6000 8000 
Current density/(Am°) 


10000 


Fig. 3. Modeling and experimental polarization curves for various H2/H20/N2 
fuel mixtures (a) 850°C and (b) 800°C. Constant H20 concentration of 4.2%. 


noted that experiments with large amounts of N3 in the fuel are 
just carried out for model validation in this study, since such 
conditions will not be encountered in actual SOFC operation. 
Fig. 3 shows a reasonable agreement between experimental and 
simulated polarization curves at 850°C, whereas at 800°C the 
model results deviates from the experimental ones at high N2 
contents (above 50% N2). 

Both modeling and experimental results indicate that the open 
circuit voltage decrease as H2 concentration decreases, which is 
more pronounced at 800°C than at 850°C. It is known that, for 
constant H2/H20 ratio, the open circuit voltage is not affected 
by the presence of inert gas [6]. The decrease in open circuit 
voltage seen in Fig. 3 is actually due to the change in H2/H20 
ratio because the water content was kept constant at 4.2%. Com- 
parison between the modeling results and experimental data in 
Fig. 3 indicates that the effects on cell performance by vary- 
ing the H2—-N2 composition becomes more significant at 800 °C 
than at 850°C, especially at lower H2 content. It is also clear 
that a limiting current density exists when N2 content reaches 
76.6% at 800°C. At these conditions, although the simulated 
polarization curve deviates from the experimental one at high 
current densities, it is worth noticing that the simulated limiting 
current density (6470 A m7?) matches well the limiting current 
density found experimentally (6430 A m7”). 

Fig. 3 highlights the effect of gas composition on cell perfor- 
mance. The gas composition affects essentially the open circuit 
voltage and anodic concentration overpotential. Fig. 4 shows 
the anodic concentration overpotential at two temperatures and 
two feed gas compositions. In this paper, the anodic concen- 
tration overpotential is not obtained from experiments. As a 
mathematical trick, if the source terms of the anode and cath- 
ode mass balance equations are set to zero, the fuel and oxidant 
concentrations are uniform throughout the electrodes and equal 
to the concentrations at the electrode/gas chamber interfaces. 
In this way, the effects of concentration overpotential could be 
separated. Since N3 is inert with regard to charge transfer pro- 
cesses and has little effects on the adsorption equilibrium and 
surface diffusion processes [7], the important characteristic of 
the feed gas is the H2/H20 ratio. In Fig. 4 the two feed con- 
centrations correspond to H2/H20O ratios of 18.2 (76.6% H2) 
and 4.8 (19.2%). It can be seen from Fig. 4 that at 850°C the 
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Fig. 4. Simulated anodic concentration overpotentials. 


Y. Shi et al. / Journal of Power Sources 172 (2007) 246-252 


anodic concentration overpotential decreases only slightly as the 
H2/H20 ratio increases from 4.8 to 18.2. However, at 800 °C, the 
anodic concentration overpotential shows a dramatic increase 
when the H2/H20 ratio is decreased. This very large increase 
in anodic concentration overpotential is the reason for the rapid 
deterioration of the cell voltage observed in Fig. 3(b) at 800°C 
and 19.2% H2. At lower temperatures (e.g. 800°C), the anodic 
concentration overpotential becomes more important due to a 
lower diffusion coefficient (diffusion through the porous anode) 
which is proportional to T!5. Also very important is the effect 
of multicomponent gas transport simulated here with the dusty 
gas model. The transport resistance increase when the H2/H2O 
ratio is decreased. 


3. Simulation results and discussions 
3.1. Hz and H20 concentration distributions 


Fig. 5 shows the H2 and H20 concentrations at the anode for 
the case of H2/H20 96%/4%, oxidant air, 800°C and cell volt- 
age 0.7 V. It can be seen that the concentration of Hz decreases 
in the electrode from the fuel chamber/electrode interface to 
electrode/electrolyte interface, while the opposite is true for 
H20. This is due to the H2 electrochemical oxidation at the 
reaction active sites resulting in the production of H20. The non- 
uniformity of Hz and H20 concentration distributions in Fig. 5 
is due to the asymmetric electrode used in our button-cell. It can 
also be seen that the profiles of species mole fraction distribu- 
tion vary faster near the electrolyte/anode active layer interface 
because of increased H? consumption and H20 production rates. 


3.2. Current density distribution 


Fig. 6 gives the electronic current density distributions within 
the anode and cathode along the cell thickness direction (r=0) 
for H2/H20 96%/4%, oxidant air and 800°C for various cell 
voltages. As expected, the cell voltage has a significant effect on 
the current density, which decreases as the cell voltage increases. 
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Fig. 6. Electronic current density distribution along the cell thickness direction 
(r=0) in (a) cathode and (b) anode. 


The electronic current density distribution can characterize the 
zone where the electrochemical reactions take place. Fig. 6(a) 
shows that the electrochemical reactions occurred within the 
whole cathode. It could then be deduced that the cell perfor- 
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Fig. 5. Surface plots of H? and H20 mole fractions. 
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Fig. 7. Effects of cathode area on cell performance in the base case. 


mance could be improved by adopting a cathode thicker than 
the 15 um layer used in this study for a thicker cathode would 
supply more TPB active sites. Fig. 6(b) indicates that in the case 
of the anode most of the electrochemical reaction occurs mainly 
within a 20 wm distance from the electrolyte/anode active layer 
interface. It is also seen that the electrochemical reactions occur 
not only within the anode active layer, but also to some extent 
within the anode support layer. From the difference in curve 
slopes between the anode active and support layers, one can see 
that the electrochemical reaction rate at the anode active layer 
is larger than at the anode support layer. This is due to the dif- 
ference in ionic conductivity and in reaction active area per unit 
volume. Thus, for this particular cell, it is likely that the cell 
performance can be improved by using a thicker anode active 
layer. However, since the mean pore size of the anode active 
layer is smaller than that of the anode support layer, increasing 
the active layer thickness will also increase diffusion resistance. 
Therefore, there should be an appropriate thickness for the active 
layer, which should not only depend on the active layer reactivity, 
but also on the layer microstructure. The effects of microstruc- 
ture parameters and electrode thickness on cell performance are 
discussed in the next sections. 


3.3. Effects of cathode area 


Fig. 7 shows the simulated effects of cathode area on cell 
polarization curves for cells with cathode areas of 0.9, 1.1, 1.3, 
1.54 (cathode area in button cell experiments), 3 and 5.3 cm? at 
800°C while keeping the anode area constant at 5.3 cm?. The 
average current density was determined on the basis of the cath- 
ode area. It is seen (in the zoom window of Fig. 7) that the 
cell performance is slightly better when decreasing the cathode 
area, especially for the smallest area considered here. However, 
Fig. 7 also shows that the effects of cathode area is really small 
when the cathode area varies from 5.3 to 0.8 cm”, which suggests 
that the effects of cathode area using an asymmetric electrode 
could almost be neglected in this case. However, the effect of 
cell performance due to asymmetric electrodes lies principally 
in the differences in anodic concentration polarization [5]. Fig. 8 
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Fig. 8. Effects of cathode area on anodic concentration overpotential. 


shows the effects of cathode area on anodic concentration over- 
potential with different fuel compositions at 800°C. This figure 
indicates that the anodic concentration overpotential increases 
only slightly as the cathode area increases from 0.8 to 1.54 cm? 
when the H2 content in fuel mixtures is high (e.g. 76.6% H2). 
However, when the H? content in fuel mixtures is low (e.g. 19.2% 
H2) the effect of cathode area is more pronounced and could not 
be neglected. In this situation, the anodic concentration overpo- 
tential shows a rapid increase when the cathode area is increased, 
especially for current densities above 4000 A m7”. This conclu- 
sion has also been experimentally supported in published data 


[8]. 
3.4. Effects of electrode porosity 


An appropriate microstructure is a compromise between 
allowing sufficient paths for gas transport and keeping enough 
three-phase boundaries surface areas for the electrochemi- 
cal reactions to take place. Porosity is one of the electrode 
microstructure parameters that usually affects the most the cell 
performance [4]. Usually, the ratio of porosity to tortuosity is 
used as a lumped parameter for representing the transport resis- 
tance. Thus, here the effects of porosity can represent the effects 
of this lumped parameter since the tortuosity was kept constant 
in the simulation. 

Fig. 9 shows the effects of cathode and anode (both support 
and active layers) porosities on cell performance when operat- 
ing at 0.7 V. Note that in this figure only one porosity is changed 
at a time. That is, the curve showing the effect of anode poros- 
ity, was determined a constant cathode porosity (36%) and vice 
versa for the effect of cathode porosity. The values of 36% for 
the porosity were determined from mercury porosimetry [1]. The 
temperature was set at 800°C and all other parameters were the 
same as those in the base case. From Fig. 9, it is seen that the 
cell performance is slightly better when decreasing the cath- 
ode porosity or increasing the anode porosity. According to the 
calculation method to determine the TPB area per unit volume 
(Stpp ), increasing the porosity leads to a decrease in the Srpg, but 
reduces mass transport resistance by increasing both the Knud- 
sen and molecular effective diffusivities [9]. For instance, our 
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Fig. 9. Effects of electrode porosity on cell performance. 


calculations have shown that when the porosity is varied from 0.2 
to 0.5, the TPB area per unit volume of the electrode is reduced 
by a factor of 4, while the molecular diffusivity doubles and the 
Knudsen diffusivity increases about four times. Thus, for a thick 
electrode, as the anode in the present case, where concentration 
polarization may becomes limiting at a higher current density, 
increasing the porosity appropriately would improve the cell 
performance. On the other hand, for a thin electrode (e.g. cath- 
ode here) where concentration polarization is not a major issue, 
increasing porosity would make the cell performance worse. 


3.5. Effects of mean particle size, mean pore size and 
electrode thickness 


Fig. 10 shows the combined effects of mean particle size 
and electrode thickness when the cell voltage is kept at 0.7 V 
and the cell temperature kept at 800°C. The changes in cell 
performance when varying the cathode thickness, anode active 
layer thickness and anode support layer thickness are shown in 
Fig. 10(a)-(c), respectively. Note that for each of these three 
figures, the layers whose thickness are not varied (e.g. anode 
active and support layers in Fig. 10(a)), are assumed to have 
a constant thickness equals to: 680 um for the anode support 
layer, 15 um for the anode active layer and cathode layer, and 
20 um for the thin-film electrolyte layer. In addition, the layers 
whose mean particle radius is not varied (e.g. particle radius 
of anode active and support layers in Fig. 10(a)) are assumed 
to have a constant value equals to: 0.6 um for the anode sup- 
port layer and cathode layer, 0.4 um for the anode active layer. 
When the mean particle size of the electrode decreases, the TPB 
area per unit volume increases and the activation polarization 
is then expected to decrease. However, it should be noted that 
when the electrode mean particle radius decreases, the electrode 
mean pore radius also decreases, therefore increasing the con- 
centration polarization because it becomes more difficult for the 
reactant gas to diffuse through the electrode due to the reduction 
of binary diffusion and Knudsen diffusion coefficients. 

Fig. 10(a) shows the combined effects of the cathode mean 
particle size and thickness on cell performance. It is seen that 
the cell performance could be improved by adopting a thicker 
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cathode and a smaller mean particles radius, and this result 
was also supported by experimental data [10]. However, if the 
particle size is very small (e.g. <0.2 wm), increasing the cath- 
ode thickness will lead to an increase in cathodic concentration 
polarization. As a result there would be no benefit to increase 
the cathode thickness above around 30 um. At a fixed cathode 
thickness, the cell performance is better with smaller mean par- 
ticle size, because of the increase in TPB area per unit volume. 
From the previous discussion, the cell used in this work could 
be improved by increasing its cathode thickness and reducing 
its corresponding particle size. 

Fig. 10(b) shows the combined effects of the anode active 
layer mean particle size and thickness on cell performance. It is 
seen from Fig. 10(b) that the cell performance can be optimized 
with smaller particle size and an appropriate active layer thick- 
ness. It should be noted that the effects of active layer thickness 
is much larger for mean particle radii larger than 0.6 um. Also, 
as the mean particle size reduces, Fig. 10(b) shows that there 
is an appropriate active anode layer thickness yielding a maxi- 
mum current density. For example, for a mean particle radius of 
0.2 um, the appropriate active anode layer thickness is around 
60 pm. This highlights the opposite effects that the mean particle 
size can have on the concentration and activation overpotentials. 

Fig. 10(c) shows that, as far as cell performance is concerned, 
an anode support layer thickness of 680 um, as adopted in the 
base case, is not necessary. The cell performance would be very 
similar with a thickness of 300 um. As for the cathode and 
the anode active layer, reducing the particle size has usually 
a beneficial effect on the cell performance. Yet, for very small 
particle sizes, concentration polarization may become dominant 
and an anode support layer too large would not be good. We 
should, however, keep in mind that the main function of the 
support anode layer is to ensure the mechanical strength of the 
PEN and thus its thickness should be kept at a relatively large 
value. To optimize the cell performance, the particle size of the 
anode support layer cannot be too small to prevent gas diffusion 
limitation. 

The above discussion emphasized the importance of opti- 
mizing the cell microstructure parameters to achieve the highest 
performance. However, the electrode thickness and particle 
size are both limited by considerations such as manufacturing 
techniques or layer mechanical strength. The present model- 
ing results should therefore, only be used as a basic guide for 
electrode design. 


4. Conclusions 


The two-dimensional isothermal mechanistic PEN model of 
an anode-supported solid oxide fuel cell for mixed gas, which 
was developed, calibrated and validated in part I of the present 
study, was further validated experimentally for extended oper- 
ating conditions by changing temperatures, anode and cathode 
gas compositions. The validation results showed that the model 


agrees well with the experimental results for a wide range of 
operating conditions. 

As expected both the modeling and experimental results show 
that the cell performance is better when operating at higher tem- 
peratures, using pure oxygen as oxidant and using a fuel with 
high hydrogen content. 

Simulation of the electronic current density distribution gives 
an indication of the extent of the zone where the electrochemical 
reaction takes place. The results with our button cell suggest that 
the electrochemical reaction occurred within the whole cathode, 
the whole anode active layer, and also to some distance within 
anode support layer. Thus, it could be deduced that the cell 
performance could have been improved by adopting a thicker 
cathode or a thicker anode active layer. A sensitivity analysis 
on the anode and cathode porosities was then performed. This 
analysis indicated that for a thick electrode, as the anode in the 
present work, increasing the porosity would improve the cell per- 
formance, while for a thin electrode, as the cathode in current 
case, increasing porosity would have a negative effect on cell per- 
formance. The effects of electrode microstructure were further 
analyzed by investigating the combined effects of the thickness 
and particle size of each electrode layer (cathode, anode active 
layer and anode support layer). The modeling results indicated 
that the cell performance for our button cell could be optimized 
by adopting a thicker cathode with small particle size, a thicker 
anode active layer (but not too thick) with smaller particle size, 
and a thinner support layer. However, the electrode thickness and 
particle size are both limited by factors such as manufacturing 
technique or mechanical strength. Nonetheless, the model pre- 
sented here can describe the effect of electrode microstructure 
on cell performance and can thus be a useful tool for electrode 
design. 
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